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It is often assumed that for treating numerical (or experimental) data on continuous transitions 
the formal analysis derived from the renormalization group theory can only be applied over a narrow 
temperature range, the "critical region"; outside this region correction terms proliferate rendering 
attempts to apply the formalism hopeless. This pessimistic conclusion follows largely from a choice 
of scaling variables and scaling expressions which is traditional but very inefficient for data covering 
wide temperature ranges. An alternative "extended scaling" approach can be made where the 
choice of scaling variables and scaling expressions is rationalized in the light of well established high 
temperature series expansion developments. We present the extended scaling approach in detail, 
and outline the numerical technique used to study the three-dimensional 3d Ising model. After a 
discussion of the exact expressions for the historic Id Ising spin chain model as an illustration, an 
exhaustive analysis of high quality numerical data on the canonical simple cubic lattice 3d Ising 
model is given. It is shown that in both models, with appropriate scaling variables and scaling 
expressions (in which leading correction terms are taken into account where necessary), critical 
behavior extends from T c up to infinite temperature. 
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INTRODUCTION 

Understanding the universal critical behavior observed 
at and near continuous transitions is one of the major 
achievements of statistical physics; the subject has been 
studied in depth for many years. It is generally con- 
sidered however that the formalism based on the elegant 
renormalization group theory (RGT) can only be applied 
over a narrow temperature range, the "critical region", 
while outside this region correction terms proliferate so 
attempts to extend the analysis become pointless. In 
fact this pessimistic conclusion follows largely because 
the traditional choices of scaling variables and scaling 
expressions are poorly adapted to the study of wide tem- 
perature ranges. 

The expressions for critical divergencies of observables 
Q(T) near a critical temperature T c and in the thermo- 
dynamic (infinite size) limit are conventionally written 



Q(T) = C Q t-«(l + F Q (t)) 
with the scaling variable t defined as 
t = (T— T c )/T c 



(1) 



(2) 



and where Fq(t) represents an infinite set of confluent 
and analytic correction terms [l[ 



F Q {t) =a Q t e + b Q t + 



(3) 



The exponents q, the confluent correction exponent 9, 
and many critical parameters such as amplitude ratios 
and finite size scaling functions, are universal, i.e. they 
are identical for all members of a universality class of sys- 
tems. When the RGT formalism is outlined in textbooks 



or in authoritative reviews such as those of Privman, Ho- 
henberg and Aharony Q or Pelissetto and Vicari Q the 
scaling variable is defined as t from the outset. However, 
because t — > oo at infinite temperature, when t is chosen 
as the scaling variable the correction terms in Fq (t) each 
individually diverge as temperature is increased. It in- 
deed becomes extremely awkward to use the expressions 
in Eq.Q outside a narrow "critical" temperature region. 
A "critical-to-classical crossover" has been invoked (e.g. 
Refs. |4, 5]) with the effective exponent 7 c ff (/?) tending to 
the mean field values as the high temperature Gaussian 
fixed point is approached. The crossover appears as a 
consequence of the definition of the exponent in terms of 
the thermodynamic susceptibility and the scaling vari- 
able t. There is no such crossover when the extended 
scaling analysis described below is used. 

Although this is rarely stated explicitly, there is noth- 
ing sacred about the scaling variable t; alternative scaling 
variables r can be legitimately chosen and indeed have 
been widely used in practice, see e.g. Refs. (r3-Hl). 

Temperature dependent prefactors can also be intro- 
duced in the scaling expressions on condition that the 
prefactor docs not have a critical temperature depen- 
dence at T r . 



An "extended scaling" approach [12Ml6l | has been in- 
troduced which consists in a simple systematic rule for 
selecting scaling variables and prefactors, inspired by 
the well established high temperature series expansion 
(HTSE) method. This approach is a rationalization 
which leads automatically to well behaved high temper- 
ature limits as well as giving the correct critical limit 
behavior. 

Here we give a general discussion of this approach. 
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We outline the relationship to the RGT scaling field 
formalism. As an illustration of the application of the 
rules, known analytic results on the historically impor- 
tant S — 1/2 Ising ferromagnet chain in dimension one 
(for which the critical temperature is of course T c = 0) 
are cited. Simple extended scaling expressions for the 
reduced susceptibility, the second moment correlation 
length, and the specific heat are exact over the entire 
temperature range from zero to infinity. An exact sus- 
ceptibility finite size scaling function is exhibited. 

The S = 1/2 nearest neighbor Ising ferromagnet on 
the simple cubic lattice is then discussed in detail. This 
model is among the principal canonical examples of a sys- 
tem having a continuous phase transition at a non-zero 
critical temperature. In contrast to the two-dimensional 
2d Ising model, in three dimensions no exact values are 
known for the critical temperature or the critical expo- 
nents. We analyze high quality large scale numerical data 
which have been obtained for sizes up to L = 256, cover- 
ing wide temperature ranges both above and below the 
critical temperature [13, The numerical technique 

is outlined. An analysis using the extended scaling ap- 
proach provides compact critical expressions with a min- 
imum of correction terms, which are accurate (if not for- 
mally exact) over the entire temperature range from T c 
to infinity and not only within a narrow critical regime. 
(The 3d Ising, XY, and Heisenberg ferromagnets have 
been discussed in Ref. fl3l]). 



DEFINITIONS OF VARIABLES 

We study the S = 1/2 nearest neighbor interaction 
ferromagnetic Ising model on the Id chain and on the 
simple cubic lattices of size L 3 with periodic boundary 
conditions. The Hamiltonian with nearest neighbor in- 
teractions of strength J is 



H 



(4) 



with the sum over nearest neighbor bonds. As usual we 
will use throughout the normalized inverse temperature 
(3 = J/kT. 

The observables we have studied are as follows: 
(i) The variance of the equilibrium sample moment, 
which is equal to the non-connected reduced susceptibil- 
ity 



X ((3,L)=N(m 2 ) = (l/N)J2(S i -Sj) 



(5) 



h3 



where m is the magnetization per spin m = (1/N) J^. Si, 
N = L d . 

(ii) The variance of the modulus of the equilibrium 
sample moment, or the "modulus susceptibility" 



Below T c , Xmod(/5, L) tends to the connected reduced sus- 
ceptibility in the thermodynamic limit, and 



(H) (/3, L) = Vx(/?,i)-Xmod(/3,L)/VA (7) 

tends to the thermodynamic limit magnetization 
(m)(/3,L) at large L. 

(iii) The specific heat which is equal to the variance of 
the energy per spin U (j3, L) 



C v (f3 7 L)=N ((U 2 )-(U) 2 ) 



(8) 



where U = (1/N) &i ' Sj with the sum over near- 
est neighbor bonds. We can note that x(P) an d C v (f3) 
have consistent statistical definitions in terms of thermal 
fluctuations. The experimentally observed susceptibility 
contains an extraneous factor f3. 

The thermodynamic limit second moment correlation 
length is defined % HI] by 

£ 2 (ft oo) =^2 (/3,oo)/2dx(/3,oo) (9) 

where the second moment of the correlation function is 

M oo) = (i/N)J2(dA -Sj) (10) 

with rij the distance between spins i and j, summing to 
infinity. 

When the "thermodynamic limit" condition L 3> 
£(/?, oo) holds all properties become independent of L and 
so are identical to the thermodynamic limit properties. 

For general L, the Privman- Fisher finite size scaling 
ansatz for an observable Q can be written 2(| 21 



Q(fi,L)/Q(p,oo)= (11) 
F Q (L/t(p,oo)) + G Q (L/Z(P,oo))) 

The functions Fq (x), Gq (x) are universal. F x (x) must 
tend to 1 when x ^> 1, and must be proportional to a; 2-1 ' 
when x <C 1. We are aware of no generally accepted 
explicit expressions for Fq (x) valid over the entire range 
of x. 



EXTENDED SCALING 

In the extended scaling approach [l2|, EH a systematic 
choice of scaling variables and scaling expression prefac- 
tors is made in the light of the HTSE. Basically, an ideal 
HTSE corresponds to the power series 



(l-y)-i = l + qy+(q(q + l)/2)y 2 + 



(12) 



Xmod(M = A((m 2 )-(H) 2 ) 



(6) 



When a real physical HTSE has the form 

Q (x) =C Q (l + a x x + a 2 x 2 H ) (13) 

with a general structure similar to but not strictly equiv- 
alent to that of Eq. (jT2"j) and a prcfactor Cq which can 
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be temperature dependent, the asymptotic limit is even- 
tually dominated by the closest singularity to the origin 
(Darboux's first theorem 22j ) leading to the critical limit 
Q{x) = Cq(1 — x)~ q . The appropriate critical scaling 
variable is 1 — x, and deviations of the series in Eq. fTB")) 
from the pure Eq. ([T2"|) form correspond to confluent and 
analytic critical correction terms. 

The extended scaling prescription consists in identify- 
ing scaling variables and prefactors such that each series 
is transposed to a form having the same structure as 
Eq. (|T3"|) , with the prefactor defined so that the first term 
of the series is equal to 1. 

The HTSE spin 5 = 1/2 expressions for the reduced 
susceptibility and the second moment of the correlation 
can be written generically in the form 0, 0, HH 



X((3) — 1 + a\x + a 2 x 2 + a 3 x 3 + 



and 



M2O8) = hx 



b 2 x 2 + b 3 x 3 



(14) 



(15) 



where a; is a normalized variable which tends to 1 as 
(3 — > p c and to zero when /3 — > 0. 

For ferromagnets, (e.g. 0, [111 [23j]) possible natural 
choices for x are x = (3/f3 c or x = tanh/3/ tanh/3 c . Scal- 
ing variables for x(/3) arc r = 1 — x = 1 — /3/j3 c or 
t = 1 — x = 1 — tanh/3/tanh/3 c . The former is stan- 
dard when T c is non-zero; when T c = 0, it is convenient 
to use x — tanh/3 (as j3 c = 00, tanh/3 c = 1). 

For /j, 2 (j3) the Eq. (fTB"]) form with the same x can be re- 
trieved by extracting a temperature dependent prefactor 
bi x so as to write 



At 2 (/3) = b lX {l + (62/61)0; + (6 3 /6i)x 2 



(16) 



The critical expressions for the reduced susceptibility 
and the second moment correlation length can then be 
written 



x(&oo)=C x t-t(1 + F x (t)) 



(17) 



(c.f. Eq.([T])) and from the relation Eq.© between \x 2 and 



ttP,^) = C i X 1 i 2 T-»{l + Ft(T)) 



(18) 



with the temperature scaling variable r = 1 — x and the 
standard definitions for the critical amplitudes C x and 
Cg. The x(/3,oo) expression has been widely used; the 
£(/?, 00) ex pres sion is specific to the extended scaling ap- 
proach fl2lll3| . The F functions contain all the confluent 
and analytic correction to scaling terms [H, l24j 



Fq(t) =a Q T e +b Q T + 



(19) 



It is important that r tends to 1 at infinite temperature 
(whereas t tends to infinity) ; the Fq (t) thus remain well 
behaved over the entire temperature range. There are 



exact closure conditions for the infinite temperature limit 
t — v 1 : C x (l + FJl)) = 1 and Q//3 C 1/2 (1 + F s (l)) = 1 
(orQ/(tanh/3 c ) 1 /2(i + F c (l)) = l). 

One can define temperature dependent effective expo- 
nents (introduced by |25j): 

7ee(r) = 51og X (^)/51ogT (20) 

see Refs. 0, For the correlation length, 

u eS (r) - d\og(m/P 1/2 )/d\ogT (21) 

is the extended scaling definition for f e g. 

For a spin S = 1/2 Ising ferromagnct on a lattice where 
each spin has z neighbors, the high temperature limit 
of the effective exponents defined by Eqns. [20] and [21] 
are 7 e ff(l) = z/3 c and ^ e ff(l) = 7cff(l)/2. A comparison 
between these values and the critical exponents 7 and 
v gives a good indication of the overall influence of the 
correction terms. If the leading confluent correction term 
in Eq. (|19[) dominates then 7 ff (1) — 7 ~ a x 6, and v c r (1) — 
v sa a(d. An analysis along these lines of j e s for Ising 
systems with large z was sketched out in Ref.[2^|. The 
case of general S is discussed in Appendix A. For all near 
neighbor Ising ferromagnets on sc or bec lattices covering 
the entire range of spin values S = 1/2 to S — 00 (which 
are all in the same 3d universality class), see Ref. 0], 
7 ff (1) and f e ff (1) differ from the critical 7 and v by a few 
percent at most. For both observables, the total sum of 
the correction terms is weak over the entire temperature 
range. 

It should be noted that traditional and widely used 
finite size scaling expressions 

QG8, L)/L^ = F Q (L^(T - T c )/T c ) (22) 

assume implicitly scaling with the scaling variable t. As a 
general rule these expressions should not be used except 
in the limit of temperatures very close to T c ; they rapidly 
becomes misleading and can suggest incorrect values of 
the exponent v if global fits are made to data covering a 
wider range of temperatures. The extended scaling FSS 
expressions [l2| 

Q{P 1 L)/(LT 1 / 2 ) q ^ = F Q ([LT 1 ' 2 ) 1 ^^ -T c )/t) 

(23) 

are valid at all temperatures above T c to within the weak 
corrections to scaling. 

For spin .5 = 1/2 Ising spins on a bipartite lattice (such 
as the Id and 3d sc lattices we will discuss below) there 
are only even terms in the HTSE for the specific heat 

SSI 

C v G8) = 1 + d r x 2 + d 2 x 4 + d 3 x 6 + ■■■ (24) 
A natural scaling expression for the specific heat is 
C v (f3) = C + C c (1 - x 2 y a (1 + F c (l - x 2 )) (25) 
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The constant term Co is present in standard analyses and 
plays an important role in 3d fcrromagncts because the 
exponent a is small. The extended scaling expression 
Eq. (f2"5"|) is not orthodox as it uses a scaling variable, r 2 = 
1 — x 2 = t(2 — r), which is not the same as the r = 1 — x 
used as scaling variable for x(P) an d £(/?)■ 



ANALYTIC RESULTS IN ONE DIMENSION 

The original Ising ferromagnet [13] consists of a system 
of S = 1/2 spins with nearest neighbor ferromagnetic in- 
teractions on a one dimensional chain. Because analytic 
results exist for many of the statistical properties of this 
system, it is often used as a "textbook" model in intro- 
ductions to critical behavior. We will use it to illustrate 
the extended scaling approach (see Ref. [HI). 

The model orders only at T = (Ref. [27|); when 
T c = the critical exponents depend on the choice of the 
scaling variable. Baxter [iij] states : " [in one dimension] 
it is more sensible to replace t = (T — T c )/T c by t = 
cxp(— 2/3)"; with this scaling variable the exponents are 
7 = 1, v = 1, a = -1 [a = -vd when T c = Ref. 

Expressions for £(/3) and x{P) m the infinite-size limit 
are readil y ca lculated following standard HTSE rules (see 
e.g. Ref. [19|). The reduced susceptibility HTSE can be 
written as 

x(P) = 1 + 2(tanh /3 + tanh 2 P + tanh 3 /? + ...) (26) 

and the HTSE for the second-moment of the correlation 
is 

fx 2 (P) = 2(tanh p + 2 2 tanh 2 P + 3 2 tanh 3 p + . . .) (27) 

The second-moment correlation length is then given by 
Eq. ([9]) with z = 2. Using the power series sums 



essentially identical for P > 2.5 but arc quite different at 
higher temperatures.) 

The internal energy per spin is just 



1709) = tanh/3 



so the specific heat 



C V (P) = cosh" 2 P 



(32) 



(33) 



Though not immediately recognizable these can all be 
re-written in precisely the form of the extended scaling 
Eqns. 1171181251 with the choice x = tanh/3 so r = (1 — 
tanh P); 

X(P) = 2(1 - tanh/3)- 1 (1 - (1/2)(1 - tanh/3)) (34) 



tanh 1/2 P 
1 — tanh P 



(35) 



and 



C V (P, oo) = (1 - tanh 2 /?) 



and so with the same critical exponents 7 



1, v 



T together with critical amplitudes C x — 2, 



(36) 

1. 

1, 

C v — 1, Co = 0. There are no analytic corrections to £(/3) 
or to C V (P) and there is only a single simple analytic 
correction to x{P)- There are no confluent corrections. 
Note again that these expressions are valid for the entire 
temperature range from T = to T = 00. 

The finite size scaling function can also be considered. 
With periodic boundary conditions the finite size reduced 
susceptibility for a Id sample of size L is 



and 



^ (l-y) 

n=l y 1,1 



VnV= yiv + 1) 



(28) 



(29) 



the exact expressions for reduced susceptibility and cor- 
relation length are thus 



and 



X(P) = cxp(2/3) 
m = i(exp(4/3)-l) 1/2 



(30) 



(31) 



(It can be noted that the "true" correlation length is 
£truc — — l/log(tanh/3). The two correlation lengths are 



X (P,L)=exp(2p) 



1 - tanh L P 



1 + tanh^ P 
The finite size scaling function is 

X (P,L)/x(P,oc) 
= tanh( J L/2^(/3, 00)) (1 + 7 _2 G X (L/2£(P, 00)) 



(37) 



(38) 



The simple principle expression 

F x {L/i{p, co)) = tanh(L/2£(/3, oo)) (39) 

is exact. 

The higher order term G x (x) in Eq. p9p is numerically 
tiny even for small L. We have not found an analytic 
expression but it can be fitted rather accurately by 

G x (x) = O.W&x 2 (l + tanh(-0.565x 1 18 )) (40) 
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HIGH DIMENSION LIMIT 



and 



For the Ising ferromagnet in the high dimension hy- 
percubic lattice limit d — > oo, with r(/3) = 1 — (3//3 c the 
reduced susceptibility and the correlation length are 



u t = t(i + Cl t + c 2 t 2 + ...) 



(46) 



and 



l/2 T -l/2 



(41) 



(42) 



exactly over the entire temperature range above T c ; the 
exponents 7 = l,v = 1/2 are of course the mean field 
exponents. These expressions again follow the extended 
scaling form given above Eqns . H~7| 1 81 including the square 
root prefactor in £(/3, 00), with no correction terms. 

In this high dimension limit the specific heat above T c 
is zero (a = 0). 

In dimensions above the upper critical dimension but 
not in the extreme high dimension limit the extended 
scaling approach has been used successfully to identify 
the main correction terms in the reduced susceptibility 

M- 

Thus analytic expressions for models both in the low 
(Id) and high (d — > 00) dimension limits follow the ex- 
tended scaling forms. This reinforces the argument that 
these forms can be considered to be generic and should 
be used at leading order also for intermediate dimensions, 
where confluent correction terms and small analytic cor- 
rection terms must be allowed for. 

In practice (e.g. Refs. 0, H E3, El 0]) analyses of 
x(/3) have long been carried out using r as the scaling 
variable rather than t. There are analogous advantages 
in scaling £(/3) with Eq. (fTS]) , which contains the generic 
{P/ (3c) 1 / 2 (or (tanh/3/ tanh^) 1 / 2 ) prefactor. We sug- 
gest that this form of scaling expression for £(/3) could 
profitably become equally standard. 



RGT FORMALISM AND EXTENDED SCALING 

In the standard RGT finite size scaling formalism 0, H| 
the free energy is written 

T% h, L) = -F sing (/3, h, L) + J- rcg (/3, h, L) (43) 

where the singular part encodes the critical behavior and 
the regular part is practically L independent. Then 

J- Sing = L- d F (u h L^ d+2 ^l 2 1 utL 1 /^ + (44) 
VuL -^) Fu (u hL ( d + 2 -^, u t lM v ) +■■■ 

with the scaling fields Uh and u t having temperature de- 
pendencies 



where h is the magnetic field. The two series in t are 
analytic. 

Ignoring for the moment the confluent correction series 
F^, for phcnomenological couplings R 

R(/3,L) = F R (u t L 1 /») (47) 
= G R (L 1 / v t(l + c 1 t + c 2 t 2 + ---fj 

and for \ 

X(P, L) = AL 2 ~ V (1 + bit + b 2 i 2 + ■■■) (48) 
= G X (l}l v t{l + ait + a 2 t 2 + ■ ■ ■ )) 

Analyses using this formalism are carried out by intro- 
ducing a series of analytic terms in powers of t, adjusting 
for each particular case the constants a n , b n and c n and 
truncating at some power of t. 

Now consider the extended scaling scheme. As a first 
step t = (T — T c )/T c is replaced in the formalism by 
t = (T — T c )/T just as for instance in [3l|]. The variable 
t is replaced by t(1 — t) everywhere. This leaves the 
generic form of the equations unchanged but modifies 
the individual factors in the series for the temperature 
dependencies of the scaling fields. 

In the extended scaling approach a second step must 
then be made due to the [fij ' Pc) 1 / 2 prefactor in 00). 
The extended scaling FSS expressions [IH 



R(/3,L)=F r [tL l ^/3 



■1/ VR -1/2 V 



(49) 



and 



X(P, L) = (L^ 2 ) 2 ^F X (50) 



can be translated into the RGT FSS formalism in terms 
of explicit built-in leading expressions for the tempera- 
ture variation of the scaling fields. The extended scaling 
expressions without correction terms are strictly equiv- 
alent to leading expressions for the scaling fields ut and 
Uh containing specific infinite analytic series of terms in 



(1 



r(l 
2v) 



r) 



-l/2v 



(51) 



%v 2 



r 2 + 



and 



u h = ha h (l + ait + a 2 i 2 ■■■) 



(45) 



u h ~ h(l - t)-( 2 -")/ 2 (52) 
^(l + ^r-^-^W-..) 
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In the extended scaling approach these leading expres- 
sions are common to all fcrromagncts. The confluent 
correction contributions will of course still exist with the 
confluent correction terms expressed using r. Finally, 
fine tuning through minor modifications of the analytic 
scaling field temperature dependence series will usually 
be necessary to obtain higher level approximations to the 
overall temperature variation of the observables. 

Not only at temperatures well above T c but already at 
criticality the extended scaling scheme can aid the data 
analysis. For instance, quite generally the critical size 
dependence of the ratio of the derivative of the suscepti- 
bility to the susceptibility is of the form 31 , 32j 



(0X08, L)/dp)/ X {P, L) = K X L 1 ^ (1 + c u L~" + ---)+K 2 

(53) 

An explicit leading order value of the L-independent term 
can be derived from the leading order extended scaling 
FSS Eq.dEB: 



K 2 = -(2- V )/2p c 



(54) 



in a ferromagnet. This value will be slightly modified by 
a correction to scaling term. 

The extended scaling scheme can thus be translated 
unambiguously into the standard RGT FSS formalism. It 
can be considered as providing an a priori rationalization 
giving explicit leading analytical temperature dependen- 
cies of the scaling fields. At this level the extended scal- 
ing scheme provides compact baseline expressions which 
cover the entire temperature region from T c to infinity, 
accurate to within confluent correction terms and resid- 
ual model dependent analytic correction terms. 



NUMERICAL METHODS APPLIED TO THE 
SIMPLE CUBIC MODEL 

The equilibrium distributions of the parameters energy 
p(U) for finite size samples from L = 4 up to L = 256 
(16, 777, 216 spins) were estimated using a density of 
states function method. When studying a statistical me- 
chanical model complete information can in principle be 
obtained through the density of states function. From 
complete knowledge of the density of states one can im- 
mediately work with the microcanonical (fixed energy) 
ensemble and of course also compute the partition func- 
tion and through it have access to the canonical (fixed 
temperature) ensemble as well. The main problem here 
is that computing the exact density of states for systems 
of even modest size is a very hard numerical task. How- 
ever, several sampling schemes have been given for ob- 
taining approximate density of states, of which the best 
known are the Wang -Landau [33[ and Wang-Swendsen 



all the information needed. Using them one can find the 
density of states in an energy interval around the critical 
region and that is all that is required for most investiga- 
tions of the critical properties of the model. 

For the present analysis a density of states function 
technique based upon the same method as in [l8| was 
used though with considerable numerical improvements 
for all L studied here (adequate improvements to the 
L = 512 data set would unfortunately have been too 
time-consuming). The microcanonical (energy depen- 
dent) data were collected as described in [l7j]. We use 
standard Metropolis single spin-flip updates, sweeping 
through the lattice system in a type-writer order. Mea- 
surements take place when the expected number of spin- 
flips is at least the number of sites. For high tempera- 
tures this usually means two sweeps between measure- 
ments and three or four sweeps for the lower tempera- 
tures we used. Note that in the immediate vicinity of (3 C 
the spin-flip probability is very close to 50% for the 3d 
simple cubic lattice. 

We report here data on the 3d simple cubic 5=1/2 
Ising model with periodic boundary conditions. For 
L = 256, the largest lattice studied here, we have now 
amassed between 500 and 3500 measurements on an in- 
terval of some 450000 energy levels, where most sam- 
plings are near the critical energy U c . For L = 128 we 
have between 5000 and 50000 measurements on some 
150000 energy levels. For L < 64 the number of sam- 
plings are of course vastly bigger. 

Our measurements at each individual energy level in- 
clude local energy statistics and magnetization moments. 
The microcanonical data were then converted into canon- 
ical (temperature dependent) data according to the tech- 
nique in (35j | . This gave us energy distributions from 
which we obtain energy cumulants (e.g. the specific heat) 
and together with the fixed-energy magnetization mo- 
ments we obtain magnetization cumulants (e.g. the sus- 
ceptibility). 

Typically around 200 different temperatures were cho- 
sen to compute these quantities, with a somewhat higher 
concentration near (3 C particularly for the larger L so that 
one may use standard interpolation techniques on the 
data to obtain intermediate temperatures. 

Below T c the variance of the distribution of m in zero 
field, Eq.([5|). represents the non-connected susceptibility; 
the physical susceptibility in the thermodynamic limit is 
the connected susceptibility 



Xconn((3,L) = N((m 2 )-((m)) 2 



(55) 



341 methods. In 17! the various methods are described 



along with an improved histogram scheme. For work in 
the microcanonical ensemble the sampling methods give 



For finite L the distribution of m below T c is bimodal 
but always symmetrical so in zero applied field (m) = 0, 
which would suggest that supplementary measurements 
are needed using small applied fields in order to estimate 
Xconn- However under the condition L >> £ C onn(/3, 00), 
where £ corm (/3, 00) is the second moment correlation 
length below T c , the two peaks in the distribution of 
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m become very well separated and the variance of the 
distribution of the absolute value \m\ can be taken as 
essentially equal to the connected susceptibility, 



XconnCA £) = Xmod(P,L) 



(56) 



The explicit expression for £ Cfmn is complicated, see [371 ] . 
but the onset of thermodynamic limit conditions can 
judged by inspection of the finite size X mo d(A L) data. 
To estimate the ordering temperature T c we have used 
the size dependence of Ui{P,L) the kurtosis of the dis- 
tribution of p{m), frequently expressed in terms of the 
Binder parameter q(0, L). 

We have introduced [38[ an alternative parameter with 
the same formal properties as g(0, L) which involves 
Xmod- The normalized parameter W(/3,L) is defined by 

W(fi, L) = 1 — ^( X mod(/3, L)/ X (I3, L))/(tt - 2) (57) 



1.58- 
1.56- 
1.54- 
1.52- 
1.50- 
1.48- 
1.46- 
1.44- 
1.42- 
1.40- 
1.38- 



Wtf, L) = (tt «|m|) 2 /<m 2 » " 2) / (tt - 2) (58) 

The normalization has been chosen such that, as for the 
Binder parameter, W = in the high temperature Gaus- 
sian limit and W = 1 in the low temperature ferromag- 
netic limit. As W(f3, L) is also a parameter characteristic 
of the shape of the distribution p(m), it can be consid- 
ered to be another " phcnomcnological coupling" . It turns 
out that at least for the 3d Ising model the corrections 
to scaling for W(/3,L) are much weaker than those for 
g(j3,L), allowing accurate estimates of T c and v from 
scaling at criticality. The values estimated for the criti- 
cal parameters f3 c and v are in good agreement with the 
most accurate values from RGT, HTSE, and Monte Carlo 
methods S. 



3D ISING FERROMAGNET SUSCEPTIBILITY 
AND CORRELATION LENGTH 

The Ising ferromagnet in dimension three is a canonical 
example of a system having a continuous phase transition 
at a non-zero critical temperature. In 3d there are no 
observablcs which diverge logarithmically in contrast to 
the 2d and 4d models. Though there are no exact results 
for this universality class, rather precise estimates of the 
critical exponents (and the critical temperatures) have 
been obtained and improved over the years thanks to ex- 
tensive analytical, HTSE, and FSS Monte Carlo studies. 
The essential aim has been to determine as accurately as 
possible the universal critical parameters. 

Consider first the finite size scaling results at and very 
close to the critical temperature. The numerical work 
[3sj ] provided an estimate /3 C = 0.2216541(2) from inter- 
sections between curves for phenomenological couplings 
at different sizes L, using data on the Binder cumu- 
lant g(L) and on the phenomenological coupling W(L) 



FIG. 1: (Color online) Finite size corrections at the critical 
temperature. x(A=, L)/L 2 ~ v against L at f3 c adopting r\ = 
0.0368. The large black points are measured; the small red 
points are the fit, Eq. (|60[) 



[381 ]. This value is consistent with the Monte Carlo es- 
timates /3 C = 0.22165452(8) and /3 C = 0.22165463(8) 
[HI, the HTSE estimate /3 C = 0.221655(2) 0, and 
p c = 0.2216546(3) [3]. 

At criticality, the standard FSS expression for 
x(Pc,L) is 



X {Pc,L)=C x L 2 



(l + d x L- u + b[L~ 



a 2 L~ 



(59) 

For the 3d Ising ferromagnet, 8 — 0.504(8) or cj = Bjv = 
0.800(13) HH so 2uj = 1.60(3). The subleading irrelevant 
exponent is u>2 = 1.67(11) [42j so the UJ2 and 2uj terms can 
be treated together as a single effective term 6 2 L -165 . In 
what follows we will assume for convenience 8 = 0.50. 

Fig. 1 shows xific, L) / L 2 ~ v against L adopting 77 = 
0.0368 [Io| ; the finite size scaling corrections in the 
present data can be fitted by 

x(Pc, L) = 1.557L 1 9632 (1 - 0.218L" 82 - 0.256L" 165 ) 

ra (60) 

The analysis is consistent with that of [10(, 

Xific, L) = L 2 ~ v (1.559(16) - 0.37(5)L-°' 8 ) (61) 

Because of the introduction of a next to leading term, 
the fit extends to lower L. 

Fig. 2 shows partial data for the ratio 



x(L) = [(dx(p,L)/d/3)/xW,L)]p c 



(62) 



against L. On this scale the data can be well represented 
by x(L) = ~{2~i 1 )/2f3 c + K 1 L- 1 /» with fi c = 0.2216549, 
rj = 0.0368(2), v = 0.6302(1), and K t a constant, see the 



CO. 
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FIG. 2: (Color online) The normalized derivative of the sus- 
ceptibility [(dxifi,L)/dP))/x(P,L)]p c against L 1,u . The ex- 
tended scaling value for the intercept is —(2 — rj) /2/3 c to lead- 
ing order (red arrow) 



extended scaling expression Eg. ([5!?]) . This form of plot 
provides an independent estimate for v consistent with 
the values given in [9I. [ToL 



To obtain an accurate value 
for v it is important to include the non-zero intercept. 

Combining v and 77 estimates from FSS at critical- 
ity, the present data are almost consistent with the MC 
and HTSE estimates 7 = (2 - t])v = 1.2372(4) 
and 7 = 1.2371(1) Q. Both of these are from meta- 
analyses on many systems in the same universality class, 
the latter relying principally on bcc data. A recent very 
precise study of the 3d Ising universality class [39| gave 
v = 0.63002(10) and rj = 0.03627(10) so 7 = 1.2372(3) 
together with u = 9/v = 0.832(6) so 6 = 0.524(4). 

Leaving the pure FSS regime, now consider the over- 
all temperature and size dependence of x(/3, L). Assum- 
ing /3 C known, the critical exponent 7 C can be estimated 
directly and independently from an extrapolation to 
t = of the derivative 7 e g(r, 00) = d\ogx{P,L)/d\ogT 
in the thermodynamic limit conditions i.e. down to 
L-dependent crossover temperatures above which the 
%(/3, L) are independent of L. The crossover occurs when 
L w 6£(/3,oo), below which the correlation length is 
no longer negligible compared to the sample size. (As 
T — > T c below this crossover, x(/3, L) then tends to a 
constant for each L). 

There is obviously no "critical-to-classical crossover" 
as a function of temperature. The crossover would ap- 
pear automatically if the effective exponent were defined 
(e.g. Ref. jij) in terms of the thermodynamic suscepti- 
bility 



FIG. 3: (Color online) An overall plot of the effective expo- 
nent 7 c ff(r, L) fixing /3 C = 0.221655 and 9 = 0.50 for sizes 
L — 64, 32, 16 from top to bottom (black, blue, green). The 
thermodynamic limit envelope curve is clearly seen. The red 
line corresponds to an HTSE data analysis [s| [4^, in full 
agreement with the present results over the entire tempera- 
ture range except for a marginal difference near /3 C . The red 
arrow indicates the consensus value for 7(/3 c ). 



and the traditional scaling variable t through 
7th,eff = <91ogxth(/3)/dlogi 



(64) 



XthCS) = [dm{P,h)/dh] h ^ = fyffl 



(63) 



because at high temperatures Xth(P) — > P and t — > T. 

The present data for L = 64, L — 32 and L = 16 are 
of very high statistical accuracy. Again assuming p c = 
0.221655, 7eff(r, L) values in the thermodynamic limit 
conditions (which are in excellent agreement with HTSE 
data for 7 e ff(T, 00) d, S^]), can be extrapolated satisfac- 
torily to r = assuming 7 e ff (t, 00) = 7 C + a\T + ■ ■ ■ , 
Fig 3. The fit provides an estimate 7 = 1.239(1), almost 
compatible with the HTSE @ and FSS [l(| estimates. 

The fluctuations in the plot for L = 64 in Fig. 4 are an 
indication of how sensitive these plots are to the slightest 
noise in the original data. The temperature region in the 
far right of Fig. 4 for L = 64 corresponds to a region 
of energy levels measured at least 500,000 times. At the 
other end the energy levels were measured more than 
1,000,000 times. Data for still higher L arc not shown 
as the fluctuations become more marked; unfortunately 
these higher L data cannot be used to refine the estimate 
of 7. The 7 estimate with the present method is sensitive 
to the value assumed for 9. The 7 estimate would become 
incompatible with the consensus value if one assumed 
significantly higher values for 8, such as 0.54 (estimates 
of 9 are reviewed in 

An advantage of this 7 e ff(i", 00) technique is that it is 
free from the problem of finite size corrections to scaling, 
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FIG. 4: (Color online) As Fig 3, blowup of the small r region. 




FIG. 5: (Color online) The normalized susceptibility 
x(/?,£)t 7 against r" assuming 7 = 1.239 and 6 = 0.50. Sizes 
L = 256,128,64,32,16,8 from top to bottom (black, red, 
green, blue, olive, orange). The excellent fit (yellow) to the 
thermodynamic limit envelope data corresponds to Eq. (f65l . 



although the Wegner thermal corrections to scaling must 
be taken into account as above. It can be noted also that 
this is a direct measurement of 7 rather than an indirect 
estimate through a combination of v c and 2— r\ c estimates 
as is the case for FSS. 

Fig. 5 shows the data for L = 16 to L = 256 in the 
form of a normalized plot, xiPi L)t 7 against t 50 as- 
suming 7 = 1.239. Again it can be seen by inspection 
at which point for each L the curves leave the thermody- 
namic limit envelope curve which is L independent. With 



FIG. 6: (Color online) The normalized correlation length 
£(/3, oo)t i/ /3 1/ ' 2 against t 9 in the thermodynamic limit assum- 
ing v = 0.630 and 9 = 0.50. Raw HTSE data provided by P. 
Butera HEl]. 

the scaling expression Eq. (|49|) and using the data at the 
various L but only in the thermodynamic limit, the fit 

X{/3, oo)r 7 = 1.106 (1 - 0.080r e - 0.016t) (65) 

gives the values of the critical amplitude, C x — 1.106(5), 
and the coefficient of the leading conformal and analytic 
correction terms, a x = —0.080(3) and b x = —0.016(3), 
read directly off the plot in Fig. 5. These values are fully 
consistent with but more precise than earlier estimates 
from HTSE, C x = 1.11(1) and a x = -0.10(3) Q, see 
[l5j ]. It can be seen that the extended scaling expression 
with only two leading Wegner correction terms gives a 
very accurate fit to the data over the whole temperature 
range above the critical temperature. 

If exactly the same data were expressed using t = (T — 
T c )/T c as the scaling variable rather than r, because r = 
i/(l + t) one would have to write 

X (/3,oo)= (66) 
1.106t _1 - 239 (l + 1.239t + 0.1466t 2 - 0.0373< 3 + ■ • • 
-0.080i ' 5 + 0.0495i 15 - 0.037K 2 5 + • • • 
-0.016t + 0.016t 2 - 0.016< 3 + ■ ■ • ) 

Remembering that t diverges at infinite T, each of the 
correction terms in the sums is individually diverging at 
high temperatures. Manifestly it is considerably more 
efficient to scale x(/3, 00) with r rather than with t. 

We have made no correlation length measurements. 
However we have carried out an extended scaling 
parametrization of HTSE thermodynamic limit second 
moment correlation length £(/3, 00) data supplied by P. 
Butera 
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FIG. 7: (Color online) The extended scaling effective expo- FIG. 8: (Color online) The leading order extended scaling 
nent v{r) against r e in the thermodynamic limit assu ming plot x(L,T)/(L(T/T c ) 1/2 ) 2 ' v against ((LT 1/2 ) 1/v t) 
and 6 = 0.50. Raw HTSE data provided by P. Butera @,[42] V ' 



Fig. 6 shows a plot of the normalized correlation length 
£(/?, co)t u ft 1 / 2 against t 8 assuming v = 0.630 and 6 = 
0.50. The data can be fitted well by the extended scaling 
Wegner expression with two leading terms only 

£{fi, oo)t v /3^ 2 = 1.074/T 1 / 2 (1 - 0.120r°- 5 + 0.051r) 

(67) 

1/2 

(note that here the critical amplitude is C^//3 C )■ The 
same equation provides the temperature dependence of 
the effective exponent defined by 



^eff 



OS.oo) = 91og(e(/?,oo)//3 1 / 2 )/aiogr (68) 



sec Fig. 7. The effective exponent varies only by a few 
percent over the whole range from T = T c to T — oo. ft 
is clear that the /3 1 / 2 prefactor is an essential part of the 
temperature dependence of the correlation length. The 
compact relation Eq. (|6"T|) is very useful as it allows finite 
size scaling analyses of the entire data set for %(/3,L). 



FINITE SIZE SCALING 

The extrapolation in Fig. 5 concerns only data in the 
thermodynamic limit condition for each L. With Eq. (j6"?l) 
in hand we can plot all the data and not just the points in 
the thermodynamic limit condition by appealing to the 
Privman-Fisher relation [45| , Eq. (fl"2j) . 

As a first step we ignore corrections to scaling and 
draw, Fig. 8, the leading order extended scaling FSS 12 1 
plot for the susceptibility 



X (L,T)/(LT^ 2 ) 2 ^ = F x ((LTV2)i/v) 



On the scale of the plot the scaling is already reasonable 
for all T above T c . 

The conformal correction can then be introduced : 

x(fi,L)/x(fi,oo) (70) 
= F x (L/£G0, co)) (1 + a x L~"G x (L/ffl, oo))) 

The function F(x) must have limits F{x) — > 1 at large 
x and F(x) ~ x 2 ^ 11 for small x. An explicit compact 
ansatz which gives these limits automatically is 

F x {x) = ((l-eM-bx^h)" (71) 



where x — L/£(J3, oo). In the critical limit x <C 1, 
F x (x)=b a (L/aP,™)) 2 - v - 



(72) 



By convention G x (0) = 1. Fig. 9 uses the tempera- 
ture dependence of the thermodynamic limit correlation 
length, Eq. (|67|) . and the thermodynamic limit suscepti- 
bility, Eq. (|6"5")) . to scale the data for all L and all (3 using 
Eq.CLED for x{P,L). 

The principle scaling function F(x) and the leading 
correction scaling function G{x) were extracted from the 
data. With the numerical constant 2 — 77 fixed at 1.963, an 
accurate effective functional form for the principal scaling 
function is 

F x (x) = [1 - cxp(-0.4179a; 1 - 963 / 1 - 262 )] 1 - 262 (73) 

On the scale of the figure F(x) with these fit values (a = 
1.262, & = 0.4179) is indistinguishable from the overall 
curve in Fig. 9. By comparing data at small L with data 
at large L the correction to scaling function can also be 



(69) 



estimated. A fit gives a x w —0.22 and 

G x {x) w cxp (-0.038a; 2 - 5 ) 



(74) 
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FIG. 9: (Color online) The Privman-Fisher scaling plot 
x(P,L)/xW,oo) against L/£_(P,oo). L = 256,128,64,32,16, 
(black, red, green, blue, cyan) 




inf 



the scaling form Eq. (|72|) has in principle only one free pa- 
rameter, a. Remarkably, when the other parameters are 
known, the FSS crossover function can be encapsulated 
in one single parameter. 

The overall scaling function expression covers all L and 
all T above T c . The principle scaling function Eq. ([7T|) 
contains only one free parameter; it resembles the finite 
size scaling form which has been used for the 2d Villain 
model |l4[ . Previous expressions for principle finite scal- 
ing functions @, HtJ , in particular for the 3d Ising model 
[H, were in the form of infinite series in exp(— a;) and so 
contained many fit parameters. 

It would be of interest to study other members of the 
same family of models in order to see if the compact form 
of scaling function Eq. (|7"Tj) is generally valid, and how the 
universality is expressed in the parameters a and b. 

Even below T c it has been noted that there should be 
a relationship between the non-connected reduced sus- 
ceptibility and the non-connected correlation length . 
The extended scaling gives explicit leading order predic- 
tions for the asymptotic relations both above and below 
T c between the finite size non-connected reduced sus- 
ceptibility x(/3, L) and the finite size non-connected cor- 
relation length £(/3,L). As we have seen, in the limit 
£(/3,L)/L«l 



xtfMLT 1 ' 2 )*^ ~ (£(/?, L)/L) 



2-r, 



(76) 



while in the opposite limit £(/3,L)/L ^> 1 the predicted 
relation is 

X (/3, L)/(LT^ 2 ) 2 -^ ~ L )/ L ){i/d)(d-2+ v ) (7?) 



For the case of the 2d square lattice Isin 
confirm both these relationships 
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model the data 
. Unfortunately, 
as we have no data here for the finite size £(/3, L) either 
above or below T c we cannot check the relationship. 



SUSCEPTIBILITY ABOVE AND BELOW T c 



FIG. 10: (Color online)The leading correction scaling function 
G(L/£(P, oo)). Black squares: measured; red circles: fit 



Fig. 10 shows the correction scaling function G(x) to- 
gether with the ad hoc Gaussian fit. 

These FSS functions are universal to within metric con- 
stants 0. 

In the same critical limit, from the definitions above 
X(/9, oo) = C x t~ 7 and oo) = Qr _ " so with 
X(fic, L) = C' I?~ n in the large L limit, 



b a = C x C 2 f n /C x 



(75) 



The ratios of susceptibility amplitudes and of leading 
correction factors above and below T c are universal. The 
standard reduced susceptibility for the region above T c 
has been discussed; for T above and below T c we will 
plot the modulus susceptibility Eq. (J6j) multiplied by |r| 7 
as a function of \t\ 6 with exponent values fixed at 7 = 
1.239,6* = 0.50, Fig. 11. 

By definition Xmod becomes equal to the connected 
reduced susceptibility below T c in the thermodynamic 
large L limit. Extrapolating the data corresponding to 
this limit to |t| =0 we find to leading order 



-1.239 



1 



10.50 



(78) 



The amplitudes C % , C x and are known from critical 
and thermodynamic limit measurements respectively, so 



with C Xi _ 



0.241(2) and a Xi _ = -0.82(5). Tak- 
ing into account the normalization factor for Xmod, the 
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FIG. 11: (Color online) The normalized modulus susceptibil- 
ity Xmod('r)|r|" > ' as a function of |r| . The upper set of curves 
corresponds to T > T c and the lower set to T < T c . In both 
cases the sizes are L = 64, 32, 16 (black, red, blue). 

present estimates for the amplitude ratio and the cor- 
rection amplitude ratio arc C Xi +/C Xj _ = 4.67(3) and 
a x . + /a x ,_ = 0.111(10). The amplitude ratio is consistent 
with previous Monte-Carlo estimates, 4.75(3), 4.72(11) 
and 4.713(7), OH 0, EH- The present correction am- 
plitude ratio estimate is however significantly lower than 
a field theory value 0.315(13) @,Hlj|. 

SPECIFIC HEAT ABOVE AND BELOW T c 

The specific heat is intrinsically difficult to analyze be- 
cause of the strong regular term Co and the small value 
of the critical exponent a (see Eq. ([23|) ). It turns out in 
addition that there are strong and peculiar finite size cor- 
rections. On the other hand the statistical precision of 
the specific heat data is very high; data for L = 512 were 
included in this analysis. The general leading form of 
the envelope data in the thermodynamic limit condition 
is assumed to be (C v (j3,L) - C ) \r 2 \ a = C c (l + a c |r 2 | e ) 
where r 2 = 1 — {(3/j3 c ) 2 . The amplitudes are C C ^ + ,C C - 
and a Ci +,a c ._ above and below T c respectively. Here a 
is fixed at 0.110, which is the expected value from the 
relation a = vd — 2 with v = 0.630. The regular term Co 
is assumed to be temperature independent; the estimate 
Co = —30.9 is obtained from the overall fit discussed be- 
low. It should be underlined that the extended scaling 
variable is r 2 not r. 

In the high temperature range (down to f3 w 0.2) the 
data can be compared to data points derived by dir ectly 
summing the HTSE terms up to n = 46 from Rcf. [52j. 
Point by point agreement is better than to 1 part in 10 3 . 



As a first step we plot the raw log C^t^L) above T c 
against logT2, Fig. 12. The thermodynamic limit data 
for different L can be clearly observed but the points fall 
on a curve rather on a straight line even down to very 
small r 2 ; this is because no Co term has been allowed 
for. Next, we plot log(C(r 2 ,£) — Co) against logr 2 , as 
in Fig. 13 for various trial values of Co- In Fig. 13 
with Co = —30.9 the envelope data now lie on a straight 
line of slope —0.110 for the lower range of t 2 (and the 
larger L). We make a Privman-Fisher finite size scal- 
ing plot of (C v (t 2 ,L) + 30.9) / (Cfo, oo) + 30.9) against 
L/£,(t2,oo) with (C„(t 2 , oo) + 30.9) taken from the ex- 
trapolated envelope for small r 2 and the measured enve- 
lope curve for higher r 2 , fitted to an explicit function for 
(C„(t 2 , oo) + 30.9), Fig. 14. The thermodynamic limit 
correlation length is taken from Eq. fTJTl) . 

In the finite size limited L < £(t 2 ,oo) region the nor- 
malized specific heat shows a strong peak, in contrast to 
the regular FSS crossover observed for the susceptibility. 
The quality of the global fit is sensitive to the value cho- 
sen for the regular term Co as the correct choice for this 
parameter is essential to obtain an i-independent peak 
height in Fig. 13. Once Co is fixed fine adjustments are 
made to the correction terms so as to obtain an L- and 
T-independent flat plateau in the left hand side thermo- 
dynamic limit region. 

An excellent global Eq. (fT2| FSS fit is obtained taking 

a(r 2 ,oo)= (79) 
-30.9 + 29.85T 2 - - n (1 + 0.12t 2 5 + 0.014t 2 ) 

The optimal value Co = —30.9(5) can be compared 
with previous estimates : —33.3(24) [53j and —27.85(80) 
0. 

The normalized (C 1 ,(r 2 , oo)+30.9)r 2 11 is shown in Fig. 
15 where the nearly linear thermodynamic limit envelope 
is obvious. 

The C„(/3, oo) from Eq. fgtJl) together with the peaked 
FSS curve (for which we have no explicit algebraic ex- 
pression) provide an accurate representation of the spe- 
cific heat at all temperatures above T c and for all sizes 
L. This is in contrast to previous analyses of MC data 
which were made in terms of truncated series of terms. 

The ratios of critical amplitudes and of leading cor- 
rection amplitudes above and below T c arc universal. 
The data show critical amplitudes C c .+ = 29.9(1) and 
C c _ = 55.4(2) above and below T c , Fig. 16. With 
the extended scaling definition, C c . + = A c 2 a /f3 2 where 
A c is the amplitude using the standard definition. The 
present C c ^ result is in very good agreement with the 
HTSE estimate A c = 1.34(1) given in Ref. [§] which cor- 
responds to C Ci+ = 29.4(3). The present estimate for 
the amplitude ratio (which is definition independent) is 
C Cj +/C c ._ = 0.540(4), consistent with e-expansion and 
field theory values of 0.524(10) and 0.541(14) respectively 
0, and with the most recent MC values 0.532(7) and 
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FIG. 12: (Color online) Raw log C„(/3, L) against logT2 
log(l - (/3//3 c ) 2 ). The sizes are from left to right L 
512, 256, 128, 64, 32, (black, red, green, blue, cyan) 



FIG. 14: (Color online) The specific heat finite size scal- 
ing fit. The ratio (C„(/3, L) - Co) / (C„(/3, 00) - Co) against 
£(/3,oo)/L with Co = -30.9. Sizes are L = 512,256, 128,64, 
(black, red, green, blue) 
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FIG. 13: (Color online) log(C„(/3,L) + 30.9) against logr 2 . 
The sizes are from left to right L = 512, 256, 128, 64, 32, 16, 
(black, red, green, blue, cyan, magenta). The dashed line has 
the slope -0.110. 

0.536(2) [iil, For the correction amplitudes the data 
indicate (Fig. 15 and 16) a c . + ~ 0.12 and a c ,_ ss —0.23, 
so a Ci+ /a c ._ ~ —0.52 and a c , + /a Xi + ~ 1.4. These values 
can be compared with field theory estimates, 0.96(25) 

and 0.95(10) respectively @,H3,|E3|. ( Jt should be noted 
that the a c values in our notation correspond to a c a in 
the notation of Refs. [HoL l5ll|.) 

We cannot carry out a full FSS analysis below T c as 
we lack information on the correlation length. 
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FIG. 15: (Color online) (C V {/3,L) - C )r£ against r 2 = 1 - 
{(3/ '/3 C ) 2 with a — 0.110 and Co = —30.9 for all temperatures 
T > T c . Sizes L = 256, 128, 64, 32, (black, red, green, blue). 



CONCLUSION 

We have applied the extended scaling approach to the 
analysis of two canonical Ising fcrromagnct models : the 
historic S — 1/2 fcrromagnet on a Id chain, and the 
5 = 1/2 fcrromagnet on the simple cubic lattice. 

For the Id model, with the scaling variables r = 
(1 — tanh/3) for the susceptibility and the correlation 
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FIG. 16: (Color online) (C v (/3, L) - C ) \T 2 \ a against |r 2 | = 
|1 - (/3//?c) 2 | with a = 0.110 and Co = -30.9. The lower set 
of curves corresponds to T > T c and the upper set to T < T c . 
Sizes L = 256, 128, 64, 32, (black, red, green, blue). 

length and T2 = (1 — tanh/3) 2 , all the analytic thermo- 
dynamic limit expressions are of precisely the extended 
scaling form over the entire temperature range from zero 
to infinity, with no confluent corrections, Eqns. I34[ 1351 

M 

An appropriate scaling variable for reduced suscepti- 
bility and second moment correlation length in a ferro- 
magnetic Ising model with a non-zero ordering tempera- 
ture is t = 1 - f5/p c = (T — T c )/T, not the traditional 
t = (T — T c )/T c . An exhaustive analysis of high qual- 
ity numerical data for the 3d Ising model demonstrates 
that the reduced susceptibility and the second moment 
correlation length can be represented satisfactorily over 
the entire temperature range above T c by compact ex- 
pressions containing two leading Wegner correction terms 
only : 

x(/3, oo) = 1.106r -1 - 239 (1 - 0.080t - 5 - 0.016r) (80) 
and 

£(£, oo) = 1.074/3 1 / 2 r-°- 630 (1 - 0.109r - 5 + 0.039r) 

(81) 

For the specific heat on a bipartite lattice (such as the 
sc lattice) the appropriate extended scaling variable is 
t-2 = 1 — (f3//3 c ) 2 . The data from T c to infinite tempera- 
ture can be fitted accurately by 

C v (fi,oo)= (82) 
-30.9 + 29.85t 2 -°- 110 (l + 0.12r 2 a5 + 0.014t 2 ) 

We give explicit finite size susceptibility scaling func- 
tions for the two models. The principle Id susceptibility 



scaling function 

X (P, L)/ X (P, oo) = tanh(L/2£(/3, oo)) (83) 

is exact. The principle 3d susceptibility scaling ansatz 

x(P,L)/ X (P,oo) = [l-exp(-fe(L/C(/3,oo))( 2 -")/ a )] Q 

(84) 

with a — 1.262, b — 0.4179 fits the data to high precision. 
This form where two parameters encapsulate the finite 
size scaling crossover from the region L ^> £(/3, oo) to the 
region L <C £(/3, oo) might well be of generic application. 

The critical parameters can be estimated by combining 
the data in the thermodynamic limit L ^> £oo(/3) with 
the data in the finite size scaling region L <C £oo(P)- 
The results provide complementary estimates for critical 
amplitudes and critical amplitude ratios. 

The aim of this work is however not so much to im- 
prove on the already very accurate existing estimates for 
universal critical parameters in the intensively studied 
ferromagnetic 3d Ising model, but to explain the ratio- 
nale leading to an optimized choice of scaling variables 
and scaling expressions for covering the whole temper- 
ature range up to infinite temperature. Here we spell 
out in detail for two canonical examples, the Id and 3d 
Ising ferromagnets, an "extended scaling" methodology 
for studying numerical data taken over the entire temper- 
ature range without restricting the analysis to a narrow 
"critical" temperature region near T c . Scaling variables 
and scaling expressions are chosen following a simple un- 
ambiguous prescription inspired by the well established 
HTSE approach. Using these and allowing for small 
leading Wegner correction terms where necessary, criti- 
cal scaling expressions for x(/3, oo), £(/3, oo) and C v (j3, oo) 
remain valid to high precision from T c right up to infi- 
nite temperature. Residual analytic correction terms are 
either strictly zero (in Id) or very weak (in 3d). 

The approach can readily be generalized to other less 
well understood systems. 

APPENDIX A : GENERAL SPIN S 

Standard expressions for the reduced susceptibility and 
the correlation length for ferromagnets as defined in [t| 
are for general spin S 

x(T) = C^std)T-^(l + F x (r)) (85) 

and 

t(T)=Cl(std)T-»(l + Fs(T)) (86) 

The extended scaling prescription consists in transpos- 
ing each HTSE expression such that it takes the form 
of a series in a variable x, having leading term 1 and 
multiplied by a prefactor. In the case of a finite critical 
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temperature Ising ferromagnet with r = 1 — /3//3 c , the 
critical amplitudes are then defined through 

X (0,oo)=C x {es)r-'r(l + F x (r)) (87) 

(c.f. Eq.flD) and 

e(Aoo) = Q(e S )^ 2 T-^(l + F £ (r)) 



For general Ising spin 5, dimension d, and a lattice 
with z nearest neighbors the extended scaling critical am- 
plitudes are 



C s x {es) 



{{S + l)/3S)C^(std) 



(89) 



and 



Cliea) 



C?(std)/ (z(l + S)/3 c /6dS) 1/2 (90) 



The definitions of the effective exponents are unal- 
tered. 

With these normalizations the physical significance of 
the critical amplitudes becomes much more transparent. 
Ref. Q lists the standard critical amplitudes as functions 
of S for sc and bec lattices. In Table AI we compare 
these values with those obtained using the above defini- 
tions. The extended scaling values are close to 1 for all S; 
the differences (Cg(es) — l) which can be read directly 
from the Table are a quantitative indication, model by 
model, of the amplitude of the ^-dependent correction 
terms within Fq. 




FIG. 17: (Color online) C^(es) plotted against C x (es), where 
C^(es) and C x {es) are spin S dependent extended scaling sus- 
ceptibility and correlation length critical amplitudes. Black 
points sc lattice, red points bec lattice. See text, Ref. [ij, and 
Table AI. 



[9J, 0.65(5) obtained by the RG in the perturbative fixed- 
dimension approach at sixth order 55j, and 0.65 from the 
e expansion to second order (Ho| . 



TABLE I: Values of the critical amplitudes for spin S with 
the standard definitions, Eqns. 1551 and 1551 Ref. 0] compared 
with values using the extended scaling definitions Eqns. l89l 
and ESI 



1/2 



3/2 



5/2 



S 
sc 

C x {std) 1.127 0.682 0.545 0.482 0.443 0.418 0.307 

Csistd) 0.506 0.458 0.443 0.436 0.432 0.430 0.423 



C x (es) 1.127 
C«(es) 1.075 

bec 



1.023 0.981 
1.003 0.979 



0.964 0.949 
0.967 0.960 



0.941 0.922 
0.957 0.945 



C x {std) 


1.042 


0.622 


0.497 


0.438 


0.404 


0.383 


0.282 


Cz(std) 


0.469 


0.426 


0.411 


0.405 


0.401 


0.399 


0.394 


Cx(es) 


1.042 


0.933 


0.894 


0.876 


0.867 


0.861 


0.845 


C\(es) 


1.023 


0.953 


0.927 


0.915 


0.909 


0.905 


0.895 



If the corrections to scaling up to infinite tempera- 
ture are dominated by the leading (confluent) term then 
Cf(es) - 1 a x (S) and Cf (es) - 1 a ( (S). The univer- 
sal ratio a<:(S)/a x (S) « (Cf(e«)-l)/(Cf (e*)-l). From 
Fig. 13 which shows the data from the Table, we can 
estimate a^(S)/a x (S) « 0.65(2) (with a small offset cor- 
responding to the next-to-leading correction) . This com- 
pares favorably with the estimates 0.76(6) from HTSE 



APPENDIX B: ISING SPIN GLASS 

It can be noted that in the case of the Ising Spin 
Glass the energy scale of the interactions is fixed by 
(Jfj) not by (J) as in the ferromagnetic case ((J,j) is 
zero in a symmetric interaction distribution spin glass). 
From an obvious dimensional argument the normalized 
spin glass "temperature" should be T 2 /(Jfj)- It has long 
been recognized that for the spin glass the HTSE ex- 
pressions contain even terms only (i.e. an expansion 
in (/3//3 c ) 2 or (tanh/3/ tanh/3 c ) 2 rather than in /3//3 c ) so 
the appropriate scalin g va riable is T sg = 1 — (f3//3 c ) 2 or 
1 — (tanh/3/ tanh/3 c ) 2 [57H60l |. The argument presented 
above for the ferromagnet can be repeated mutatis mu- 
tandis on this basis; the extended scaling expressions for 
x(/3, L) and £(/3, L) in spin glasses are the same as those 
for the ferromagnet (Eqns. HTl and fl"8|) but with (/3//3 c ) 2 
substituted for /3/j3 c everywhere [l2| . 

Unfortunately the great majority of publications on 
spin glasses have used t as the scaling variable which is 
quite inappropriate except for a very restricted range of 
temperatures near T c . One consequence is that many 
published estimates of the exponent v in spin glasses are 
low by a factor of about 2 (see the discussion in [32| ) . 



16 



ACKNOWLEDGEMENTS 

We would like to thank Paolo Butera for generously 
providing us with tabulated data sets and for helpful 
comments. We thank V. Privman for an encouraging 
comment. This research was conducted using the re- 
sources of High Performance Computing Center North 
(HPC2N). 



[25] 



F. J. Wegner, Phys. Rev. B 5, 4529 (1972). 

V. Privman, P.C. Hohenberg and A. Aharony, in Phase 

Transitions and Critical Phenomena, edited by C. Domb 

and J.L. Lebowitz , Vol. 14 (Academic Press, New York, 

1991). 

A. Pelissetto and E. Vicari, Phys. Rep. 368, 549 (2002). 
E. Luijten, H. W. J. Blote, and K. Binder, Phys. Rev. 
Lett. 79, 561 (1997). 

Y. Garrabos and C. Bervillier, Phys. Rev. E 74, 021113 

(2006) . 

M. Fahnle and J. Souletie, J. Phys. C 17 L469 (1984). 
S. Gartenhaus and W. S. McCullough, Phys. Rev. B 38, 
11688 (1988). 

J.-K. Kim, A. J. F. de Souza and D. P. Landau, Phys. 
Rev. E 54, 2291 (1996). 

PButera and M. Comi, Phys. Rev. B, 65 144431 (2002). 
Y. Deng and H. W. J. Blote, Phys. Rev. E 68, 036125 

(2003) . 

M. Caselle, M. H asenbusch, J. Phys. A 30, 4963 (1997) 
|hepdat/970T0"07] .4.75(3) . 

I. A. Campbell, K. Hukushima, and H. Takayama, Phys. 
Rev. Lett. 97, 117202 (2006). 

I. A. Campbell, K. Hukushima, and H. Takayama, Phys. 
Rev. B 76, 134421 (2007). 

H. G. Katzgraber, I. A. Campbell and A. K. Hartmann, 
Phys. Rev. B 78, 184409 (2008). 

I. A. Campbell and P. Butera, Phys. Rev. B 78 024435, 
(2008). 

K. Hukushima, LA. Campbell and H. Takayama, Int. J. 
Mod. Phys. C 20, 1 (2009). 

R. Haggkvist, A. Rosengren, D. Andren, P. Kundrotas, 
P. H. Lundow and K. Markstrom, J. Stat. Phys. 114 455 

(2004) . 

R. Haggkvist, A. Rosengren, P. H. Lundow, K. Mark- 
strom, D. Andren and P. Kundrotas, Adv. Phys. 56 653 

(2007) . 

M. E. Fisher and R. J. Burford, Phys. Rev. 156, 
(1967). 

E. Brezin, J. Phys. (Paris) 43 15 (1982). 
P. Calabrese, V. Martin-Mayor, A. Pelissetto, and E. 
cari, Phys. Rev. E 68, 036136 (2003). 

J. G. Darboux, J. Math. Pures Appl. 4, 377 (18 78). 
P. Butera and M. Comi, |arXiv:hep-lat /0204007 

F. J. Wegner, in Phase Transitions and Critical Phenom- 
ena, vol 6, ed C Domb and M S Green (New York: Aca- 
demic Press) (1976). 

J. Kouvel and M. E. Fisher, Phys. Rev. A 136, 1626 
(1964). 



583 



Vi- 



[26] G. Orkoulas, A. Z. Panagiotopoulos, and M. E. Fisher, 
Phys. Rev. E 61, 5930 (2000). 

[27] E. Ising, Z. der Physik 31, 253 (1925). 

[28] R. J. Baxter, Exactly solved models in statistical me- 
chanics, Academic Press (1982). 

[29] G. A. Baker and J. C. Bonner, Phys. Rev. B 12, 3741 
(1975). 

[30] B. Berche, C. Chatelain, C. Dhall, R. Kenna, R. Low, 
and J. -C. Walter, J. Stat. Mech. P11010 (2008). 

[31] M. Hasenbusch, A. Pelissetto, and E. Vicari, Phys. Rev. 
B 78, 214205 (2008). 

[32] H. G. Katzgraber, M. Korner and A. P. Young, Phys. 
Rev. B 73, 224432 (2006). 

[33] F. Wang and D. P. Landau, Phys. Rev. Lett. 86 2050 

(2001) . 

[34] J. -S. Wang and R. H. Swendsen, J. Stat. Phys. 106 245 

(2002) . 

[35] P. H. Lundow and K. Markstrom, Cent. Eur. J. Phys. 7 
490 (2009). 

[36] K. Binder, Z. Phys. B: Condens. Matter 43, 119 (1981). 
[37] H. Arisue and K. Tabata, Nucl. Phys. B 435, 555 (1995). 
[38] P. H. Lundow and I. A. Campbell, Phys. Rev. B 82, 

024414 (2010). 
[39] M. Hasenbusch, Phys. Rev. B 82, 174433 (2010). 
[40] J. Salas and A. D. Sokal, J. Stat. Phys. 98,551 (2000). 
[41] R. Guida and J. Zinn- Justin, J. Phys. A: Math. Gen. 31, 

8103 (1998). 

[42] K. E. Newman and E. K. Riedel, Phys. Rev. B 30, 6615 
(1984). 

[43] P. Butera, private communication. 

[44] P. Butera and M. Comi, Phys. Rev. B 58, 11 552 (1998). 
[45] V. Privman and M. E. Fisher, Phys. Rev. B 30 322 
(1984). 

[46] M. E. Fisher, in Critical Phenomena, Proceedings of the 
51st Enrico Fermi Summer School, edited by M.S. Green 
(Academic Press, New York, 1972). 

[47] S. Caracciolo, R. G. Edwards, S. J. Ferreira, A. Pelis- 
setto, and A. D. Sokal, Phys. Rev. Lett. 74, 2969 (1995). 

[48] J. Engels, T. Scheideler, Nucl. Phys. B 539, 557 (1999). 

[49] M. Hasenbusch, Phys. Rev. B 82, 174433 (2010). 

[50] C. Bagnuls and C. Bervillier, Phys. Rev. B 24 1226 
(1981). 

[51] C. Bagnuls, C. Bervillier, D. I. Meiron and B. G. Nickel, 

Phys. Rev. B 35 3585 (1987). 
[52] H. Arisue and T. Fujiwara, Phys. Rev. E 67, 066109 

(2003) . 

[53] M. Hasenbusch and K. Pinn, J. Phys. A 31, 6157 (1998). 
[54] X. Feng and H. W. J. Blote, Phys. Rev. E 81, 031103 
(2010). 

[55] C. Bagnuls and C. Bervillier, J. Phys. A 19, L85 (1986). 
[56] M. C. Chang and J. J. Rehr, J. Phys. A 16, 3899 (1983). 
[57] R. Fisch and A. B. Harris, Phys. Rev. Lett. 38,785 
(1977). 

[58] R. R. P. Singh and S. Chakravarty, Phys. Rev. Lett. 57, 
245(1986). 

[59] L. Klein, J. Adler, A. Aharony, A. B. Harris and Y. Meir, 

Phys. Rev. B 43 11249 (1991). 
[60] D. Daboul, I. Chang, and A. Aharony, Eur. Phys. J. B 

41, 231 (2004). 



!)!) 



